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Abstract 

In this paper, we study spectral properties of generalized weighted Hilbert matrices. In particular, 
we establish results on the spectral norm, determinant, as well as various relations between the 
eigenvalues and eigenvectors of such matrices. We also study the asymptotic behaviour of the spectral 
norm of the classical Hilbert matrix. 



m 
o 

s 

in 
-i— > 



> 
in 
oo 

o 

co 
O 
m 



x 



1 Introduction 



The classical infinite Hilbert matrices 






-1 


1 

2 


1 

3 


1 





-1 


1 
2 


1 

2 


1 





-1 


1 

3 


1 

2 


1 






and H n 



/I 

2 

J_ 

3 

1 

4 



(1) 



'•7 



have been widely studied in the mathematical literature, for a variety of good reasons (see [2] for a nice 
survey of their astonishing properties). In this paper, we present results and conjectures on spectral 
properties of these matrices and related types of matrices. We first review known results in Section [2] and 
then introduce new results in Section [3] on generalized weighted Hilbert matrices of the form 







if m = n 
if m ^ n 



Our results can be summarized as follows. Theorem[T]below states a surprising property of these matrices: 
their spectral norm depends monotonically in the absolute values of their entries, a property known a priori 
only for matrices with positive entries. A second important result (Theorem [5]) is that the determinant 
of such matrices are polynomials in the square of their entries. We prove next in Lemma [5] a key relation 
between the eigenvalues and eigenvectors of these matrices, that lead to a chain of nice consequences 
(among which Corollaries Q] and [5]) . The present work finds its roots in the seminal paper of Montgomery 
and Vaughan [7], which initiated the study of generalized Hilbert matrices. 

Notations. Let p > 1. In what follows, ||y|| p denotes the £ p -norm of the vector y € C , i.e. 




and for an S x S matrix M, \\M\\ P denotes the matrix norm induced by the above vector norm, i.e. 

\My\\ P 



\M\\ P := 



sup 

|y||p=i 



In the particular case p = 2, the following simplified notation will be adopted: 

II y II 2 = ||y|| (Euclidean norm) and ||M|| 2 = ||M|| 
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Notice in addition that when M is normal (i.e. when MM* = M*M, where M * stands for the complex- 
conjugate transpose of the matrix M), the above norm is equal to the spectral norm of M, i.e. 



||M|| = sup{|A| : A <= Spec(M)} 



2 A survey of classical results and conjectures 
2.1 Hilbert's inequalities 

The infinite-dimensional matrices presented in ([lj are two different versions of the classical Hilbert matrix. 
Notice first that T x is a Toeplitz matrix (i.e. a matrix whose entry n, to only depends on the difference 
m — n), while is a Hankel matrix (i.e. a matrix whose entry n, m only depends on the sum n + to). 

The original Hilbert inequalities state (see [U p. 212]) that for u, v e £ 2 (Z;C) (resp. in u, v e £ 2 (N;C)) 
with Hull = llvll = 1, 



^ U m (Too )m,nV n < 7T resp. I ^2 u m(-ffoo)r 



< 7T 



where it cannot be replaced by a smaller constantQ. This is saying that Too (resp. i/oo) is a bounded 
operator in £ 2 (Z;C) (resp. in £ 2 (N;C)), with norm equal to ir. 

Hardy, Littlewood and Polya obtained an explicit expression for ||i?oo||p m [5 P- 227], for all values of 
p > 1: 

ll-floollp = . , , r , P>1 
sm(7r/p) 

and Titchmarsh proved in |13| that UTooUp < oo. Also, UTooUp is clearly greater than or equal to ||ifoo||pj 
as Hoc may be seen as the lower- left corner of (up to a column permutation), but no exact value is 
known for it (except in the case where p = 2" or p = 2 n /(2™ — 1) for some integer n > 1; see [5J for a 
review of the subject). 

Consider now the corresponding finite-dimensional matrices Tr and Hr of size R x R: 

!0 if to = n 2 
1 .. . and (H R ) mn = ■ l<m,n<R 
lfTO^n m + n — 1 
m — n 

The above Hilbert inequalities imply that for every integer R > 1, both 

IITrH < Tr and \\H R \\ < Tr (2) 

Clearly also, both ||Tr|| and increase as R increases, and limfl._j.oo ||Tr|| = limfl_ s . o = t 1 - 

A question of interest is the convergence speed of ||-fffl|| and ||Tr|| towards their common limiting value 7r. 
Observing that up a column permutation, H R can be seen as the lower- left corner of T2_r+1j we see that 
\\Hr\\ f° r every integer R > 1. This hints at a slower convergence speed for the matrices H R 

than for the matrices T R . Indeed, Wilf et de Bruijn (see [TS]) have shown that 

7r ~ l|gfl|l ~ 2(log(fl))^ aSi? ^°° 
whereas there exist a, b > such that 



£<*-||T*||<^ {R >2) (3) 



1 Notice that Hilbert had proved originally these inequalities with 2ir instead of it; the optimal constant was found later 
by Schur. 
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We will prove these inequalities at the end of the present paper. The lower bound has already been proved 
by H. Montgomery (see [TC]) and it has been conjectured in pTO] and independently by Montgomery, that 
the upper bound in the previous inequality is tight, i.e. that 

„„, „ c logi? 
t- T« #— asi?^cx> 

We also provide some numerical indication of this conjecture at the end of the paper. 
2.2 Toeplitz matrices and Grenander-Szego's theorem 

We review here the theory developed by Grenander and Szego in [3] for analyzing the asymptotic spectrum 
of Toeplitz matrices. In particular, we cite below their result on the convergence speed of the spectral 
norm of such matrices. 

Let (c r , r€ Z) be a sequence of complex numbers such that 

^|c r |<oo (4) 

and let us define the corresponding function, or symbol: 

f{x) = c r exp(zra;), x € [0, 2ir] 

Because of assumption made on the Fourier coefficients c r , / is a continuous function such that /(0) = 
/(27r) (equivalently, / can be viewed as a continuous 27r-periodic function on K). 

Let now Cr be the Rx R matrix defined a 

(Cfl)m,n = Cm-n, 1 <m,n < R 

The following fact can be verified by a direct computation: for any vector u E C R such that ||u|| 2 = 

J2l<n<R \ U n\ 2 = 1) 

<-27T 

u*C R u = / /(a:) |0(a;)| 2 dx (5) 
Jo 

where f(x) is the above defined function and <j>(x) — — ]= X)i<n<i?, M « exp(z(n — l)x). 



Let us now assume that Cr is a normal matrix (i.e. CrC* r — C* r Cr)] this is the case e.g. when / 
is a real-valued function (in which case Cr is Hermitian, i.e. C* R = Cr). As ||u|| = 1, we also have 
Jq V |^>(a;)| 2 dx = 1, which implies that 

\\C R \\< sup \f(x)\=:M 

for any integer R > 1. Grenander and Szego proved in (31 p. 72] the following refined statement of 
the convergence speed of the spectral norm. If / is twice continuously differentiable, admits a unique 
maximum in xq and is such that f"(xo) ^ 0, then 



M-\\C R \\~f(x )-f(x + ^) 

as R — > oo. 



^ 2 |/"K)I 
2R 2 



The above theorem does however not apply to Hilbert matrices of the form T Rl as the harmonic series 
J2 r >i h diverges, so condition ^ is not satisfied. Correspondingly, the symbol associated to these 
matrices is the function 

„. . — exp(irx) + exp(— irx) „ v-^ sm(rx) , . . ^ . 

/(a:) = 5^ ^ '- ^ ^ = _2z^^^ = i (x-7r), a;e]0,27r[ 
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while by Dirichlet's theorem, /(0) = f(2n) = 0; / is therefore discontinuous, but relation still holds 
in this case and allows to deduce Hilbert's inequality: 

\\Tr\\ < sup |/(a;)| = vr 

xG[0,2Tr] 

However, relation ([5]) alone does not allow to conclude on the convergence speed towards ir. 

In general, the problem of evaluating the convergence speed of the spectral norm is a difficult one when 
/ attains its maximum at a point of discontinuity. An interesting matrix entering into this category 
has been the object of a detailed study by Slepian in (see also [U] for a recent exposition of the 
problercH). It is the so-called prolate matrix, defined as 



(PR)m,n = Prn-n, l<m,n<R, where p r = 



' sin(27rmr) [f r ^ Q 

2irw if r = 



where < w < | is a fixed parameter. Here again, we see that condition (U) is not satisfied. The symbol 
associated to this matrix is the function 

/ i ■ \ n sin(27rwr) 
Jw(x) = }_^Pr exp(irx) = 2nw + 2^ cos(rx) = tt I[o,2™]u[27r(x-u0,2,r](a0 

rGZ r>l 

for all x G [0, 27r]\{27T'u;, 27r(l — w)}. In this case, we again have for any integer R > 1 

||.Pr|| < sup |/„(o;)|=7r and lim H.F'rII = 7r 

a:e[0,27r] R ^°° 

It is moreover shown in [T^] that for all < co < | , there exist c w ,d w > (where both c w and d^, are 
given explicitly in [T3]) such that 

tt - H-PrII ~ c w VR exp(-du,i?) 

We see here that even though the function f w is discontinuous, the convergence speed is exponential, as 
opposed to polynomial in the case of a smooth symbol. Of course, the situation here is quite particular, 
as the function f w has a plateau at its maximum value, which is not the case for the Hilbert matrix Tfj. 

2.3 Generalized weighted Hilbert matrices 

Let x = (xi, . . . , xr) be a vector of distinct real numbers, c = (ci, . . . , cr) be another vector of real 
numbers, and let us define the matrices A(x) and P(x, c) of size R x R as 

!0 if m = n 
1 -r -l ( 6 ) 
ltmf n 
Xn 

and 

{0 if m = n 
c m c n (7) 
lfmf n 

If there is no risk of confusion, the matrices A(x) and i?(x, c) will be denoted as A and B, respectively. 



2 We would like to thank Ben Adcock for pointing out this interesting reference to us. 
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2.4 A result and a conjecture by Montgomery and Vaughan 



In order to motivate the study of generalized Hilbert matrices, let us mention here both a result and a 
conjecture made by Montgomery and Vaughan. The next section will be devoted to applications of these 
questions. 

In [7], Montgomery and Vaughan showed the following result. If ^4(x) is the R x R matrix with entries 
as in ([6]), then 

P(x)|| < ~ 5 

where 8 = mfi< TOl „<.R, m ^„ \x n - x m \. 

They further showed that if -B(x, c) is the R x R matrix with entries as in ([7|) with c n = & n d 
5 n = mini< m < fl . m ^„ \x m - x n \, then 

||B(x,c)||< T (8) 

They also conjectured that the tightest upper bound is actually ||S(x, c)|| < tt. In [9J, Montgomery and 
Vaughan's result was improved to ||S(x, c)|| < 4^, but the conjecture remains open so far. 



2.5 Applications 

2.5.1 Large sieve inequalities 

Let Xi, . . . ,xr be real numbers which are distinct modulo 1. Let also ||t|| be the distance from the real 
number t to the closest integer and let 

8:= min ||a: r — x s \\ and S r :— min — a; s || 

For (a„) M+i<n<M+N an arbitrary sequence of complex numbers, we write 

S(x) := a n exp(27rma;) 

M+l<n<M+N 

A large sieve inequality has the generic form 

\S(x r )\ 2 <A(N,S) J2 l a «! 2 

l<r<_R M+l<n<M+N 

Using Hilbert's inequality @, one can show that the previous inequality holds with A (N, 8) = N+5^ 1 — 1. 
Equivalently, this says that 

if B := {exp(27rma; T .)} M+1 < n < M+JVil < r .< i j, then ||B|| 2 < A(N,S) 

Besides, generalized Hilbert inequalities of the type © are particularly useful when studying irregularly 
spaced x r (such as Farey sequences). These generalized inequalities allow to prove the following refined 
large sieve inequality: 

E ( N+ l S r 1 ) W,)| 2 < E 
l<r<_R ^ ' M+l<n<M+N 

This last result is useful for arithmetic applications, as it allows e.g. to show (see [3) that tt(M + N) — 
7r(M) < 27r(7V) (tt(N) being the number of primes smaller than or equal to N), whereas the inequality 
7r(M + N) — tt(M) < tt(N) stands as a conjecture so far. 

Another important application of large sieve inequalities is the Bombieri- Vinogradov theorem (see for 
instance pQ), which is related to various conjectures on the distribution of primes. 
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2.5.2 Other Hilbert inequalities 

In [8], Montgomery and Vaughan study variants of Hilbert 's inequality (with for instance - \ replaced 
by csc(x r — x s )), which allow them to show the following result: if X) n >i n l a «| 2 < °°j then 



f I Va„n- lt \t = V|a„| 2 (T + 0(n)) 
Jo ' 



n>l n>l 

The key idea behind the proof of the main result in [5] is the following identity: 

csc(x fc - xi) csc(xi - x m ) = csc(x k - x m ) (cot(x k - x{] + cot(a; i - x m )) 

which is of the same type as our relation (fTU)) below. Building on this, a further generalization of Hilbert's 
inequalities has been performed in where the following functional equations are solved: 

= *(x) -*{y)+ t (X ~ V l and 1 = ~ ^ + t{x)t( V ) (with t(0) = 0) 



6(x)6(y) y ' vy/ 0(x-y) 6(x)6(y) 8{x - y) 



3 New results 



3.1 Spectral norm of 5(x, c) 

We establish below a monotonicity result regarding the spectral norm of matrices of the type -B(x, c). 
Theorem 1. 7/x, x',c, c' are vectors of real numbers such that 

I u m,n 

(x, c)| < |6 m ,„(x',c )| for 1 < to,7i < R 

then 

||B(x,c)||<||B( X / ,c')|| (9) 

Remark 1. Notice for matrices Y, Z with positive entries, it holds that if < y m ,n < z rn . n for all m, n, 
then ||y|| < \\Z\\. Consider indeed the normalized eigenvector u corresponding to the largest eigenvalue 
of Y*Y: as Y*Y has positive entries, u is also positive, so ||Y|| = \\Yu\\ < \\Zu\\ < \\Z\\. The above result 
states that a similar result holds for matrices of the form -B(x, c), even though these do not have positive 
entries. 



The remainder of the present section is devoted to the proof of Theorem [TJ which we decompose into a 
sequence of lemmas. 

First of all, let us observe that the numbers a m ,„ = l/(x m — x n ) satisfy 

&k,l ^l,m — ^k,m 

(ak,i + ai !m ) for k, I, to distinct (10) 
This relation will be of primary importance in the sequel. 

Lemma 1. If k is a positive integer and 1 < 77 < R, then, denoting as B- n the matrix B with n th row 
and column removed, we obtain 

S:= Klb m ,n{B k _ n \ m = V (11) 

l<l,m<R 



Proof. Using (fTU)) . we obtain 



5' 



E 

l<l,m<R 



l.m 



0>n,l 



){B k - n )i 



Si + S2 



l<;,m<_R 
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with 

l<l,m<R l<l,m<R l<m<R 

and 

^2= cic m (? n a m ^ia n ,i{B k L n )i. m = c 2 n a nt i {B k i r ^)i,i = — c 2 , ai i?l (B k j n 1 )i,j = — Si 

l<l,m<R 1<1<R 1<1<R 

as A is antisymmetric. □ 

Lemma 2. Lei 1 < n < i? and k > 2 be an integer. Then 

(B k ) n ,n= 51 51 bn > 1 ( B -n)l,m bm - n ( Bk T ~ 2 )n,n = ~ X! X! ^ B -n)l,l (B^^n^ 

0<r<fc-2 l</,m<i?, 0<r<fc-2 l</<fl 

Proof. Notice first that 

(B ) n ,n = ^ ' 011,71! b nit n 2 ■■■ b nk _ 2t n k l bn k _ 1 ,n 

l<ni,...,nfc_i<fi 

As = 0, we may consider m, rifc-i ^ n in the above sum. For each (n\, . . . , nfc_i), define 

s = inf {t £ {2, . . . , k} | m ^ n, . . . , n t _i ^ n,n t = n} 

(where by convention, we set rife = n in the above definition). Ordering the terms in the above sum 
according to the value of s, we obtain 



(B k ) n , n - ( B -n) ')„,„ 



X] 5Z ( B -n) ni ,n r+t h ^r+i.n {B k r 2 ) ^ 



0<r<fc— 2 ni,n,.-|.i/n 

which is the first equality in the lemma. The second one follows from (fTT|) and the fact that B is 
antisymmetric. □ 

Lemma 3. Let 1 < n < R and k > 2 be an integer. Then the following holds 

- if k is odd, then {B k ) n:Tl = 

- if k is even, then (— l)^ (B k ) n ^ n is a polynomial in (bf m , 1 < I < m < R) with positive coefficients. 

Proof. Since B is antisymmetric, the first statement is obvious. The second one follows by induction 
from Lemma [2j □ 

Proof of Theorem^ Observe that since the matrix iB is Hermitian, it has R real eigenvalues fj,i, . . . , \ib. 
corresponding to an orthonormal basis of eigenvectors, so 

\\ B \\ = max \ix r \ 

l<r<R 

and for a positive integer k 

Tr(B 2k ) = (-l)Vf 

l<r<R 

Therefore, we obtain 

= lim ((-l^Tr^ 2 *))^ 

k— >oo 

and the theorem follows from Lemma [3l □ 
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3.2 Determinant of 2?(x, c) 

The following result shows that the determinant of -B(x, c) is a polynomial in bf m . 

Theorem 2. -If R is odd, then det(£?(x, c)) = 0. 
-If R = 2T is even, then 

R T T 

det(5(x,c)) = n^ £ A £ LK^, (12) 

fc=l (mi,ni)fe-B i=l (m^n^feS i=l 

where 

E := {{m l ,n i ) 1 [ \ uj =1 {m^nj = {1,...,R} and m l < m, Vi, mi < . . . < to t } 
Let us first establish the following lemma. 

Lemma 4. Let I be an integer, with 3 < I < R. Denoting by Si the set of permutations of {1, ...,/}, we 

have 

S ■= £ 0<t(1),<t(2) a <r(2) : cr(3) ' ' ' a cr(;-l) : cr(;) <V(i),o-(i) = (13) 

o-eSi 

Proof. Let us define 

ft := £ <V(1),ct(2) a <r(2),cr(3) • • • a cr(;-l),a(l) a<7(/-i),o-(0 

and 

ft '■= £ a <r(l),cr(2) a <T(2),<x(3) 1 1 1 a <T(Z-l),<r(l) Oct(/),ct(1) 

By pXJj) . we have 5 = ft + ft. Let now t e 5; be the permutation defined by r(l) = / — 1,t(2) = 
1,t(3) =2,...,t(1-1) = 1- 2,t(1) = I. We obtain 



ft — £ «CTT(l) : ffT(2) fl ffT(2),JT(3) ' ' ' a CTT(i-l),CTT(l) QfTT (Z) , CTT ( 1 ) 

= £ fflff(J-l),<r(l) a <r(l) : cr(2) ' ' ' Ct(r(i-2),tr(i-l) a cr(i),CT(Z-l) = ~ft 



which completes the proof. □ 



Proof of Theorem^ By definition, 



det(.B) = £ ( fJ ) H V. 



,fr(n) 

o-e-Sij l<n<_R 



Now, every permutation er is a product of cycles, say F\, F%, . . . ,Fk. Let us denote by n\ , n% , . . . , rik the 
respective cardinalities of these cycles and let us set 

^(-^ ) ■ ^ ^ ^Si ,S2 a S2 ,S3 ' ' ' ^ , S rli — X,S rti ^S n . ,Si 

Si ,S2,...,S ni |{si,S2,...,S n ^ }=J*i 

In the above expression for det(_B), the contribution of the permutations having these sets as support for 
their cycles is 

fc R 
(_ 1 )n 1 +n 2 +...+»*-ft JT JJ C 2 

i— 1 r— 1 

By (|13p and the fact that the main diagonal is zero, a non-zero contribution can therefore only occur 
when all cycles are of cardinality 2, which proves the theorem. □ 



Remark 2. The above statement allows to recover part of the conclusion of Lemma [3J First notice that 
by Theorem [5] and for all J C {1, . . . , R}, det(Bj), where Bj = {b^ m )^ mG j, is also a polynomial in tif . 
Define then 

ak = e n x ' 

JC{l,...,R}ieJ 

\J\=k 

where Ai, . . . , Xr are the eigenvalues of B. Notice that 

o k = det (^) ( 14 ) 

JC{1,...,R} 

\.J\=k 

Indeed, let P be the polynomial defined as P(x) — Y[i<i<R( x ~ ^0- We observe that on one hand, the 
matrix-valued version of this polynomial is given by 

p(x)= n (x-xd) = x R + J2 xR ~ k (- 1 ) k E Ii*i = * R + E 

l<i<R k=l ,7C{1 R}i£J l<k<R. 

\J\=k 

while on the other hand, 

R R 

P{x) =JJ(s-Ai) = det(xI-B) = x R + ^x R - k {-if ^ det(Bj) 

i=l fc=l Jc{l,...,fl} 

|J|=fc 

so identifying the coefficients, we obtain equality (|I4[) . This implies that a k is also a polynomial in tif m . 
Finally, for si = X^k^h-Mj we have the following recursion, also known as Newton-Girard's formula: 

Sl = i-iy-^isi-i + i-iy-Ha! 

l<i<i-l 

For example, sq = n,si — o~±, S2 — s\o~\ — 2a2, S3 = S2 o~\ — s\02 + 3fi3, etc. We therefore find by 
induction that for all k, (—l) k Tr(B 2k ) = ( — l) k S2k is also a polynomial in bf m , but this alone does not 
guarantee the positivity of the coefficients, obtained in Lemma [3] above. 

3.3 Formulas regarding the eigenvalues and eigenvectors of A(x) and £?(x, c) 

We first state the following lemma, which has important consequences on the eigenvalues of the matrices 
^4(x) and B(x, c), as highlighted hereafter. The approach taken below generalizes the method initiated 
by Montgomery and Vaughan in [TJ. 

Lemma 5. a) Let u = (ui, . . . , ur) t be an eigenvector of A(x) for the eigenvalue i^. Then for I < n < R, 
we have 

A* 2 K| 2 = a m,n(K»| 2 +2 9?K^)) (15) 

l<m<R 

b) Let u = (iti, . . . , un) T be an eigenvector o/_B(x, c) for the eigenvalue ifi. Then for 1 < n < R, we have 

^\Unf = fl l« (^dk| 2 + 2c n C m%n^)) (16) 

l<m<_R 

Proof. Clearly, (|15|) is a particular case of (fT6|) (with all c„ = 1). In what follows, we prove (fTBf directly. 

Our starting assumption is Bu = i/zu, i.e. Yli<m<R b n ,m u m = i[iu n . Taking the modulus square on 
both sides, we obtain 

fJ 2 \u n \ 2 = Y b n,mb n ,lU m Ul 
l<l,m<R 



9 



(Notice that the sum can be taken over I ^ n, m ^ n as b n>n — 0). Therefore, 



\j? \u 



n\ 2 = Cn cic m a n , m a n ,iu m ui = c 2 n (Si + S 2 ) (17) 



l<;,m<fl 



where Si corresponds to the terms in the sum with I = m, i.e. 

•51= C mOm,J u m| 2 (18) 



l<m<R 



and 

5*2 = ^ cic m a n , m a nt iu m ui 



\<l,m<R 
l^m, l^n, m^n 



As I, m, n are all distinct in the above sum, we can use f| 10[) and the antisymmetry of A gives 

&n,m ^nj — &Z,m ^n,/ T" Om,l ^n,m 



SO 



with 



5*2= Q c m (a; !m a n ,i + a mji a„, m ) u m u; = S3 + Si (19) 



l<l,m<R 
l^m, Z^n, m^n 



S3 = ^2 c l c mai, m a n! lU m Ul = ^2 h, m an,lU m Ui = ^2 a n,lUi ^2 h,mU 

l<l,m<R l<l,m<R 1<1<R l<m<R 

As u is an eigenvector of B, it follows that 

S3 = a nt iui(ifj,ui-bi lV ,u n ) 

1<1<R. 

Likewise, noticing that u is also an eigenvector of B (with corresponding eigenvalue — i/i), we obtain 

^4 ^ ^ Q>n,m ^ ^ ^m.l 1^1 ^ ^ ^n,m ( ifAUui ^m,n ^n) 

l<m<_R. l<i<_R. \<m<R 

From (Tl9l). we deduce that 



(U m U n + U m U n ) = 2 ^ Un) 
l<m<fi l<m<fl 

Now, using this together with (TIT)) and (fT5)) , we finally obtain 

M = ^ , c n ( c m a m,nl w m| + 2 C m C n d m rl 5i(u m U„)) 

l<m<fl 



which completes the proof. □ 
One of the many consequences of Lemma [5] is the following. 

Corollary 1. If C\, . . . , cr are all non-zero, then the eigenvalues o/i?(x, c) are all distinct. 
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Proof. Indeed, if in the basis of eigenvectors of B, there were two eigenvectors corresponding to the same 
eigenvalue, then it would be possible to find a linear combination of them (which is also an eigenvector) 
such that one component (say u n ) would be equal to zero. Then by (|16p . we would have 

a m,n C n C m\ U m\ = 

l<m<R 

which is impossible, given the assumption made. □ 
A more precise version of Lemma Ob) reads as follows. 

Lemma 6. Let u = v + iw be an eigenvector of -B(x, c) (with v = 3? (u), w = 3(u) ) corresponding to 
the eigenvalue ifj,, then 

^ V n = X b n,m w m + ^ c n ^2 a n,m W m {[I V m - b„ hn W n ) (20) 
l<m<fl, l<m<fl 

Moreover, if [i 7^ 0, then ||v|| = ||w||, while if fi = 0, then det(B) = 0, so one of the eigenvectors 
corresponding to this eigenvalue is real. 

Proof. Applying the proof method of Lemma [5] gives 

v l = ( X] bn ' m Wm ) = X bn, m w m + X bn ' m bnl Wm Wl = ^ + ^ 2 

l<m<R l<m<R l<l,m<R 



with 



and 



2 2 
w 

n.m m 



Km<R 



S2 = X bn ' m bn l Wm Wl = C n X ° l ° m an - m an - 1 Wm Wl = C ™ ^ 3 + 



l<l,m<R l<l,m<R 



with again 



S3 = X Cl ° m ai > m a ™>' Wm Wl = X an > lWl X b 

l<l,m<R 1<1<R l<m<R 

= X a n,l w l(l lv l-bl,nW n ) 
1<1<R 

and likewise, 

&L = X a n,rnW m ^ b m jWi= ^ a n,m W m (^ V m - b m , n W n ) 
\<m<R 1<1<R. l<m<R 

Observing that S3 — S4, we obtain the formula (|20|) . 
Finally, we have by assumption £?(v + iw) = ifj, (v + iw), so 

Bw = /iv and Bv=-/iw 

Consequently, 

|| w|| 2 = /iw T w = (-Bvfw = (B T v) T w = v T Bw = /i ||v|| 2 
so for 7^ 0, we have ||v|| = ||w||. □ 
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Finally, let us mention the following nice formula. 

Lemma 7. Let u be an eigenvector of B corresponding to the eigenvalue fj,, Then 



|2 

l<r<R l<r<R 



Proof. Let C — diag(ci, . . . , cr) and X — diag(a;i, . . . , xr). Then 

u T {XCAC - CACX) u = u T M u 
where m r , s = c r c s for r ^ s and otherwise. Therefore, 



u T M u — \ c r u r — | c r u 



- , <v| 2 

l<r<_R. KKR 



On the other hand, 

u T (XCAC - CACX) u = u T (IB - Bl) u = u T Ii/iu-i/iu T Iu = 
as u T (-B) = u T B T = (Bu) T = (-i^uf = -ifiu T . The result follows. □ 

3.4 Back to the spectral norm 

Lemma [5] also allow us to deduce the following bounds on the spectral norm of A(x). 
Corollary 2. 

max a m,n < P( x )l| 2 < 3 max Y] a 2 (21) 

~ ~ l<n<R ~ ~ l<n<R 

Proof. The left-hand side inequality is clear, as the m th column of A is the image by A of the m th 
canonical vector. For the right-hand side inequality, we use (fT6|) . choosing n such that \u n \ 2 > |u m | 2 for 
aU 1 < m < R, and = ||A||. We therefore obtain 



^H 2 |«n| 2 = a 2 m ^ n (\u m \ 2 + 2M{u n U^)) < ^ndWml' + kml' + l^l 2 ) 

l<m<_R 

A|| 2 |U„| 2 <3 <nWn\ 2 



l<m<R l<m<R 

so 



Km<R 



□ 



3.5 The classical Hilbert matrix Tr 

First of all, notice that the upper bound in equation (|21[) allows to recover to the original upper bound 
on ||Tr||, where Tr is the Hilbert matrix defined in the introduction: 

\\T R \\ 2 < max 3 V ^<3-2V^=^ 2 

Km<R ^ (m~ n) 2 ^ n 2 

Let us now come back to the convergence speed of ||Tr|| towards 7r, already mentioned in Section[2] We 
shall now prove inequality |3]), namely the fact that there exist positive constants a and b such that 

a ,. b \os(R) 

\\ t r\\ < f (^>2) 
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The lower bound can be deduced from Lemma [5] From equation ([IT))) , we indeed see that if R = 2S + 1, 
then 

S 



\t r \\ 2 < 6 y i = tt 2 - 6 y i < tt 2 - e y — 1— = ^ 

11 ^ k 2 ^k 2 ^Lk(k+1) 



so 

_ |, T I, 6 3 

^ " *" (S + l)(7r + ||T fl ||) tt(S+1) 
which is indeed of the type ^ < tt — 

Another way to prove this lower bound is to follow the Grenander-Szego approach of Section 12.21 Let us 
first recall equation ([S]): 

2tt 

2 



u*T R u= / f(x)\(ft(x)\ 2 dx 
Jo 

where f(x) — i (x—tt) for x <E (0, 27r) and <ft(x) — Yli<n<R u « exp(i(n— 1)jc), where both J*^ ^(x)! 2 <fa = 
Hull 2 = 1. Therefore, 



7T — U*iTrU 

where 



/*2tT /•27T 

/ x|0(a:)| 2 da; so tt- 11^11= inf / x \(ft(x)\ 2 dx (22) 
Jo ^eE(R) J 



E{R) = \<j)(x) = -^= y u n exp(i(n - l)x) ueC R , y |u„| 2 = 1 



2tt i 

l<n<i?. 



l<n<7? 



It remains therefore to show that the term on the right-hand side of (f2"2l is bounded below by a term of 
order 1/ R. To this end, let us consider (ft 6 -E(-R) and c > 0: 



/ |</>(x)| 2 <fa = — y* u m u n [ cxp(i(m ~ n)x) dx 

J 27T l<tZ<R J° 



2 

C \ - , , C 



< £ y ww^l E iki <f E M 2 = f 

l<m,n<R \l<n<R J \<n<R 

where we have used Cauchy-Schwarz' inequality. Setting therefore c= |,we obtain Jq \<fi(x)\ 2 dx < \. 
This in turn implies that for all (ft £ E(R), 

2tt 



x\(ft{x)\ 2 dx> x\(ft{x)\ 2 dx > - \(ft(x)\ 2 dx> n 

Jtt/R K Jtt/R Zli 

which settles the lower bound in equation ([3]). 

In order to establish the upper bound, we need to find a function (ft € E(R) such that 

^ x\(ft(x)\ 2 dx< (23) 
R 

for some constant b > 0. This will indeed ensure the existence of a vector u, the one associated to the 
function (ft € E(R), such that |u*Tru| > tt ~ °^ R , implying the result. 

In view of equation (|2"3")h our goal in the following is to find eft £ E(R) such that for both c and e small, 



2tt 



\cft(x)\ 2 dx < e (24) 



This would indeed imply that 



27T f-C />27T 

x \cft(x)\ 2 dx < c \<ft(x)\ 2 dx + 2tt / |</)(x)| 2 da; < c + 27re (25) 
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Let M and N be positive integers such that N{M — 1) + 1 < R and let 

g{x) — f exp(ima;) 



N 

ixp(imx) ) 
Q<m<M-l 



The function </> defined as <p(x) — g(x — c/2) / J J Q 27r \g(x)\ 2 dx belongs to E{R). Our claim is that for an 



appropriate choice of M and N, 4> satisfies (|24[) with both c and e small. 
We first show the following estimate on f 2 * |g(x)| 2 dx. 
Lemma 8. 

M 2N i r 2n 

< _L / | 5 ( x )|2 dx < M 2iV-l 



iV(M - 1) + 1 - 2tt Jo 
Proof. Let -ftT be a positive integer and define the polynomial 

K 



p K {t) = { £ t m ) = W 

Q<m<M-l 0<i<KT(M-l) 

Notice that clearly, = 6 m ,_R- if I + m = K(M — 1). Moreover, 

\g(x)\ 2 = \P N (exp(ix))\ 2 = £ b LN b m , N exp(i(l - m)a 

0<l,m<N(M-l) 



SO 



fZTY 

/ |ff(a;)| 2 da; = 2tt £ = 2tt £ 6/^6 

0<i<JV(M-l) 0<i<JV(M-l) 

What remains therefore to be proven is 

M 2N 2JV _ 1 

< 0JV(M-1),2JV < ^ 



JV(M-l)-i,JV — 27t6jv(M-1) : 27V 



7V(M - 1) + 1 
Using Cauchy-Schwarz' inequality, we obtain 



(Eo 



2 



v ^0<i<AT(M-l)^,A r J Pat(1) 2 _ M 27V 

VN(M-1),2N — 2-^i \n — 



N(M-l) + l N(M-1) + 1 N(M-1) + 1 

On the other hand, P2N(t) = Pi{t) P2N-i{t), so 

°N(M-1).2N = £ h,2N-l < P2iV-l(l) < M 2JV_1 

(JV-l)(M-l)<Z<jV(M-l) 

which completes the proof. □ 



We now set out to prove (|24|) . Recall that (j>{x) = g(x—c/2) j y J Q 27r |g(x)| 2 dx. As a result of the previous 
lemma, 

2 " \x< m2j ^iV(M-l) + l 1 f 2 ". 2 JV(M-1) + 1 1 f 2 *- c/2 . , |2j 
|0(z)| 2 tfa< A/2jV ; ^/ |ff(z-c/2)| 2 cfa = M2JV ; — / \g(x)\ 2 dx 



c/2 



Notice that 

i / m2 I V- /■ a 2N fsm(Mx/2) x 2A 



Q<m<M-\ 



sin(x/2) 
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Figure 1: Rescaled gap f(R) between the spectral norm of the infinite-dimensional operator and that 
of the matrix T R , as a function of R G {1, . . . , lO'OOO}. 



so 



\g(x)\ 2 dx = 2 \g(x)\ 2 dx<2 1 

Jc/2 Jc/2 



'c/2 Jc/2 

as for < x < it, sin(a;/2) > x/tt. This implies 

r 2-K-c/2 



dx 



/2 

and correspondingly 



\g{x)\ 2 dx<2 (-) 

Jc/2 ^ xJ 



2N 



dx = 2tt 



/27T V 



2N d y 



2n 



2N - 1 



2N-1 



27r \H^ dx < N(M ' 1) + 1 



1 



Assuming R > 3 and defining M 



2R 



log R 



N :-- 



M 2N 



lo S R 



27V- 1 



2iY-l 



and c := — -ff — (where [x] denotes the integer 

is satisfied with 



part of x), we verify that M(N — 1) + 1 < R (so g E(R)) and prove below that 
£ = 0(1/ R). Indeed, as M > and iV(M - 1) + 1 < M(2N - 1), we obtain 

- 1) + 1 / Af/9 1 + N(M - 1) 

M^(2Jy-l)(c/27r)^-i = (cM/2 " } M(2jV-l) " (CM/27F) 

as 1 — 2A^ < 3 — logi?. According to ([23)1 , this finally leads to 



l ~ 2N < exp(l - 27V) < - 



x|</>(x)| 2 dx < 



7re log R 2ire 3 



R 



R 



which completes the proof of the upper bound in ([3]). As already mentioned, it has been conjectured in 
[TU] that of the two bounds in ((3]), the upper bound is tight. We provide below some numerical simulation 
that supports this fact; on Figure [I] the expression 



f(R) :=(7r-||Tfl||) 



R 



logi? 



is represented as a function of R, for values of R ranging from 1 to lO'OOO: 

Detailed facts can also be established about the eigenvectors of Tr. In order to ease the notation, suppose 
that R = 2S + 1 and that Tr is indexed from —S to S. 

Lemma 9. Let u be an eigenvector of Tr, corresponding to the eigenvalue ifi, with Uq = 1 (one can 
always multiply u by a constant in order to ensure that this is the case). Then for < n < S, we have 
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Figure 2: Amplitude {|it„|, —R < n < R} of the eigenvector corresponding to the largest eigenvalue of 
T R , with R = l'OOO. 



Proof. Define v by v n — —U— n . Then 

(T B y ) . m = y: -^—= e -s-=- e - 

* J — 777 — 77 L ' — 77? 4- T7 ^ ' T77 



m - n -m + n ^ m - n 

S<n<S —S<n<S —S<n<S 



So _ 

(7fiv)_ m = V" — = (T R u) m = (-i/j,u) m = inv- 

— ' m. — n. 



m — n 

-S<n<S 



i.e. v is an eigenvector corresponding to the eigenvalue ifj,, with Vq = 1. Thus, v = u (as the eigenspace 
corresponding to ifx is of dimension 1). □ 



We finally conjecture that the following statement holds. Let u be the eigenvector corresponding to the 
largest eigenvalue /i in absolute value. Then 

\u m \ < \u n \ V0 < m < n < S 

This conjecture is confirmed numerically; on Figure [2j we represent |tt n | as a function of n £ {— S, . . . , S}, 
for S = l'OOO. 

From the theoretical point of view, the above conjecture seems also reasonable, as (— l) fe (To fc )„ >n (see 
Lemma [2]) should decrease as n increases (in absolute value). If true, this fact would therefore hold in 
the limit k — ¥ oo, which would imply the conjecture on the eigenvector. 
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